Planetary Nebulae an Eye to the Future 
Proceedings IAU Symposium No. 283, 2011 
A.C. Editor, B.D. Editor & C.E. Editor, eds. 



©2011 International Astronomical Union 
DOI: 00.0000/X000000000000000X 



Planetary Nebulae in the Magellanic Clouds 
and Local Group Galaxies 



00 



Warren A. Reid 1 2 

department of Physics and Astronomy, Macquarie University, Sydney, NSW 2109, Australia 
2 Macquarie University Research Centre in Astronomy, Astrophysics & Astrophotonics 
email: warren.reid@mq.edu.au; war@aao.gov.au 

7— I ' 

' Abstract. The Magellanic Clouds are close enough to the Milky Way to provide an excellent 
environment in which to study extragalactic PNe. Most of these PNe are bright enough to 
£Nj . be spectroscopically observed and spatially resolved. With the latest high resolution detectors 
q ' on today's large telescopes it is even possible to directly observe a large number of central 
stars. Magellanic Cloud (MC) PNe provide several astrophysical benefits including low overall 
extinction and a good sample size covering a large range of dynamic evolutionary timescales 
while the known distances provide a direct estimation of luminosity and physical dimensions. 
on : Multi- wavelength surveys are revealing intriguing differences between MC and Galactic PNe. 

Over the past 5 years there has been a substantial increase in the number of PNe discovered 
in the Large Magellanic Cloud (LMC) in particular. Deep surveys have allowed the faint end of 
the luminosity function to be investigated, finally providing a strong clue to its overall shape. 
In so doing, the surveys are approaching completeness, estimated at ~80% in the LMC (~120 
deg 2 ) and ~65% in the Small Magellanic Cloud (SMC) (~20 deg 2 ). 
i-C . The number of galaxies comprising the Local Group (LG) and its outskirts has been growing 

Oh- steadily over the past 5 years and now numbers 48. Most of the 7 newly discovered galaxies 
are dwarf spheroidal (dSph) in structure and range from 7.6 to 755 kpc from the Milky Way. 
!— i ' Nonetheless, there are no published searches for PNe in any of these galaxies to date. Apart 
from the LMC and Milky Way, the number of PN discoveries has been very modest and only one 
additional LG galaxy has been surveyed for PNe over the past 5 years. This paper provides the 
number of Local Group PNe currently known and estimates each galaxy's total PN population. 
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1. The advantages of Local Group PNe outside the Milky Way 

Most astrophysical PN parameters, including ionized nebular mass and the brightness 
and evolutionary states of their central stars, depend on accurate distances (Ciardullo 
y—{ ' et al. 1999). This is difficult in our own Galaxy due to inherent problems with variable 
t-H . extinction and lack of central star homogeneity (Terzian, 1997). Since extra-galactic 
PNe are located outside the Milky Way they have the advantage of well established 
distances. The chemical composition of individual PNe within each galaxy can be directly 
compared in order to derive the intrinsic range in galactic chemical composition as it was 
a few Gyrs ago. Indeed, the oldest progenitors of today's observed PNe did not process 
elements such as S and Ar and to some extent also O and Ne. This makes them excellent 
probes for estimating the metallicity of the ISM at the time when the progenitors were 
agglomerating in the early life of their host galaxy. Elements produced and processed 
within the PN progenitor up until the end of the AGB including He, C, N and the s- 
process elements are returned to the ISM, enriching material for the next generation 
of stars, allowing models of stellar evolution to be tested in a galactic context. Strong 
emission lines enable the measurement of electron temperatures and densities, which, 
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along with emission line fluxes, provide estimates of central star temperatures and mass. 
Accurate radial velocities can also be measured using a number of emission lines. 

These galaxies in turn represent a diverse mix of morphological types, metallicities, 
chemical evolution, star formation history and galactic evolution. Since all of these fac- 
tors directly affect the observed PNe, including characteristics of their central stars and 
surrounding nebulae, a study of these objects becomes one of galactic archeology. Oxygen 
and neon abundances measured in PNe can be related to present day abundances in the 
ISM. By plotting the galactic location of observed PNe, models of past star formation 
and galactic evolution can be created. 

2. The advantages of Magellanic Cloud PNe 

PNe in the Magellanic Clouds are able to provide key data on the physics of stel- 
lar evolution, mass loss (Iben 1995), nucleosynthesis processes, abundance gradients, 
ISM chemical enrichment and powerfully trace star-forming history (e.g. Maciel & Costa 
2003). The well-determined 50 Kpc LMC distance (e.g. Keller & Wood 2006, Reid & 
Parker 2010a), modest, 35 degree inclination angle and disk thickness (only ^500 pc, 
van der Marel & Cioni 2001), mean that LMC PNe are effectively co-located. Since dim- 
ming of their light by intervening gas and dust is generally low and uniform (e.g. Kaler 
& Jacoby, 1990), absolute nebula luminosity and size can be estimated with minimal 
correction. 

Although there are dwarf galaxies nearer to the Milky Way, few if any PNe have been 
found in them. The Magellanic Clouds provide single environments at known distances 
with luminosities up to 100 times that of the closer galaxies. Surveys in the Magellanic 
Clouds can not only identify candidate PNe but, with the benefit of known distance, 
provide initial emission-line luminosities, local environmental conditions and spatial de- 
tails. With the combination of multi-wavelength data now available, including ground- 
based spectroscopy (eg. Reid & Parker 2010a), HST imaging and UV spectroscopy (eg. 
Stanghellini ct al 2005), mid infrared data from the Spitzer Space Telescope survey of 
the LMC: Surveying the Agents of a Galaxy's Evolution (SAGE Meixner et al (2006)), 
deep YJK S photometry (eg. Cioni et al 2011), and radio observations (eg. Filipovic et 
al 2009) the astrophysical tools are finally available to understand the effects of mass, 
temperature and metallicity during the PN formation and evolutionary cycles in each 
Magellanic Cloud. 

The direct measurement and comparison of both nebulae and central stars across a wide 
evolutionary range allows the construction of extremely accurate nebula photoionization 
models. This task is more difficult to achieve in the Milky Way or more distant galaxies. 
Such powerful data contributes to our understanding of dredge-up effects, expansion 
velocities, changes according to dynamical age, nebula chemical contribution to the ISM, 
past stellar motions including those due to tidal forces and galactic evolution. 

3. PN populations 

Stellar population synthesis theory shows that the PN population size, when scaled 
to the available completeness limits, directly correlates with the visual magnitude of the 
galaxy (Magrini et al. 2003, Corradi & Magrini 2006). This relation, however, is weakly 
dependent on the age, metallicity and the degree of star formation at the time the PN 
progenitors were forming. In addition, the PN specific luminosity rate (the proportionality 
coefficient: number of PNe L~ feo; ) or "a ratio" (Jacoby 1980) may be strongly influenced 
by the morphological type of the host galaxy (Buzzoni, A., Arnaboldi, M. 2006). 
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When the theory of simple stellar populations (SSPs) (Renzini & Buzzoni, 1986 and 
Buzzoni, 1989) is applied to the Local Group, we have a basic reference frame allowing 
the evaluation of "a" for coeval and chemically homogeneous stellar conglomerates. In 
simple terms, the expected number of PN (Npjy) for a SSP of total luminosity (L^) is: 

Npjv = fi L bo; t pn (3.1) 

where 6 is the so-called "specific evolutionary flux" (which approximates to 2 x 10 -11 
L~J o/ yr _1 ; see Buzzoni (1989)) and Tpjq is the PN (emission-detectable) lifetime in years. 

In reality, however, the stellar core mass evolution directly affects the PN lifetime since 
it depends on the central star temperature (T e /y ~ 10 5 K) during the hot post AGB 
phase (see Renzini & Buzzoni, 1986). Using models from Paczynski (1971) and Vassiliadis 
and Wood (1994), Buzzoni & Arnaboldi (2006) have expressed this fuel in terms of H- 
equivalent solar mass. With an increase in the hot post AGB liftime with increasing 
dynamical ages and metallicities, the thpagb is somewhat equivalent to Fuel/lnpAGB, 
where Ihpagb is the luminosity of the central star at the onset of the PN stage. The 
luminosity-specific PN density then becomes: 

NpjV 

a = = Qt pn = &mm{T HP AGB,Tdyn} (3.2) 

Lssp 

Despite physical age constraints that exist for the SSP, Buzzoni & Arnaboldi (2006) 
have estimated a theoretical luminosity-specific PN density based on their equation (3.2) 
where 

a = 1PN/1.5 x 1O 6 L . (3.3) 

With a small bolometric correction applied to the L v of the galaxies, it is possible 
to use a to estimate the number of PNe that may exist in each galaxy. This bolometric 
correction does not amount to much more than -0.2 Mag. but appears to bring a closer to 
several independent population estimates eg. ~1,000 PNe for the LMC (Reid & Parker, 
2006a,b), and 8,000 ± 1.500 for the Andromeda galaxy (Peimbert, 1990). In order to 
compare the number of PNe currently uncovered and the likelihood of further discoveries, 
the known number of PN for each galaxy is plotted on the same scale and shown in Figure 
[TJ The a estimate for each galaxy is shown by the open circles, which have a strong linear 
correlation. The main reason for real PNe deviating away from the same correlation, even 
at lower PN density levels, is the lack of a homogeneous completeness limit. The number 
of spectroscopically confirmed and candidate PNe currently known in the Local Group 
galaxies is shown in column 6 of Table [T] (2011). Column 5 gives the number of PNe 
known at the time of the previous PN conference (IAUS234 in 2006). It can be seen that 
the largest increase in the number of candidate PNe has occurred in the LMC and in 
our own Galaxy. Other Local Group galaxies record only moderate increases, but the 
numbers expected to be found will also be much lower (see Figure [T]) due to their lower 
overall stellar populations. Column 7 provides a reference to the 2006 census whereas 
column 7 provides a reference to the most up to date and extensive census work. 

4. Multi-wavelength PN surveys and research in the Magellanic 
Clouds over the past 5 years 

Over the past 5 years there have been 16 published papers that specifically analysed 
PNe in both the LMC and SMC, 11 of which employed space-based observations. In 
addition there were individually 8 LMC and 6 SMC PN papers published, 8 of which 
used space-based data. Overall, 63% of the PNe studies in the MCs over the past 5 years 
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Figure 1. The number of PNe currently known in each galaxy as a function of the absolute 
estimated luminosity of their host galaxies. Filled circles represent the current census in each 
galaxy. Open circles represent the possible number of PNe predicted using a, based on the theory 
of simple stellar populations. 



employed space-based observations using Spitzer Space Telescope (NIR), HST (optical 
and UV), FUSE and XMM- Newton satellite (x-ray). The past 5 years has importantly 
seen major improvements to radio telescope detectors, leading to the mapping of radio 
emission from MC PNe. These detections have already become the subject of 4 papers. 

4.1. Optical observations 

The largest and most productive survey for PNe in the LMC over the past 5 years has 
been that conducted by Reid & Parker (2006a, b). This UK Schmidt Telescope survey 
of the central 25deg 2 main bar region used the best 12 best deep Ha and 6 short red 
images, obtained over a 3 year period, which were then scanned and combined using 
SuperCOSMOS. Candidate identification was carried out by assigning a separate colour 
to each digital image and then overlaying the two. This method is somewhat different to 
the usual methods of on-band/off-band [O in] imaging and wide- field slitless spectroscopy 
such as objective prism imaging. The process is fully described in Reid & Parker (2006a). 
The resulting candidates were spectroscopically followed up using a variety of telescopes 
including 2dF on the A AT, and FLAMES on the VLT leading to new LMC kinematic 
studies (Reid & Parker 2006b), a new PN luminosity function (Reid & Parker 2010a) 
and a new improved method for determining excitation class (Reid & Parker 2010b) . 

Over 2,000 candidate objects of different types were spectroscopically confirmed, in- 
cluding PNe in various stages of evolution which were released with a quality /certainty 
flag. Continuing multi- wavelength observations (see FigureEJ have secured classifications, 
reducing the number of candidates in the lowest quality (possible PN) and uncovering 
93 new PNe in the outer LMC (see Reid & Parker 2011). 

4.2. Abundances 

Using the empirical method for abundance determination which applies ionisation cor- 
rection factors (ICFs), Leisy & Dennefeld (2006) found a slight anti-correlation between 
N and O for MC PNe. This relation which is more pronounced in the SMC is probably 
due to higher efficiency of the surface enrichment through dredge-ups which is expected 
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Table 1. Known planetary nebulae populations in the local group and outlying galaxies in the 
years 2006 and 2011 showing the improvement. Local group galaxies newly identified/discovered 
in the past 5 years are indicated by an asterisk next to the name in column 1. 



Name 


Type 


Mv 


Dist. 


PNe 


PNe 


Rcf (old) 


Rcf (new) 








[kpc] 


2006 


2011 


2006 


2011 




M31 


Sb 


-21.2 


785 


2766 


2766 


Merrett 2006 




Milky Way 


Sbc 


-20.9 




2400 


3000 


Acker et al. 1996 


Parker et al. 2006: 
















Miszalski et al. 2008 


M33 


Sc 


-18.9 


795 


152 


152 


Ciardullo et al. 2004 




LMC 


Ir 


-18.5 


50 


277 


740 


Jacoby 2006 


Reid 2006a,b, 2011 


SMC 


Ir 


-17.1 


59 


105 


139 


Jacoby et al. 2002 


Jacoby 2006 


M32 (NGC221) 


E2 


-16.5 


760 


30 


45 


Ciardullo et al. 1989 


Sarzi et al. 2011 


NGC205 


Sph 


-16.4 


760 


35 


35 


Corradi et al. 2005 




IC10 


Ir 


-16.3 


660 


16 


27 


Magrini et al. 2003 


Kniazev, et al. 2008 


NGC6822 


dlr 


-16.0 


500 


17 


26 


Lcisy et al. 2005 


HM 2 et al. 2009 


NGC185 


Sph 


-15.6 


660 


5 


5 


Corradi et al. 2005 




IC1613 


dlr 


-15.3 


725 


3 


3 


Magrini et al. 2005 




NGC147 


Sph 


-15.1 


660 


9 


9 


Corradi et al. 2005 




WLM 


dlr 


-14.4 


925 


1 


1 


Magrini et al. 2005 




Sagittarius 


dSph/E7 


-13.8 


24 


3 


4 


Zijlstra 1999 


Zijlstra et al. 2006 


Fornax (E351-G30) 


dSph 


-13.1 


138 


1 


2 


Danziger et al. 1978 


Larson 2008 


Pegasus (DDO 216) 


dlr 


-12.3 


760 


1 


1 


Jacoby et al. 1981 




Leo I (DDO 74) 


dSph 


-11.9 


250 










Andromeda I 


IDsPH 


-11.8 


810 










Andromeda II 


dSph 


-11.8 


700 










Leo A 


dlr 


-11.5 


690 


1 


1 


Magrini et al. 2003 




DD 210 


dlr 


-11.3 


1025 










Sag DIGD 


dlr 


-10.7 


1300 










Pegasus II 


dSph 


-10.6 


830 










Pisces (LGS3) 


dlr 


-10.4 


810 










Andromeda V 


dSph 


-10.2 


810 










Andromeda III 


dSph 


-10.2 


760 










Leo II (Leo B) 


dSph 


-10.1 


210 










Cetus* 


dSph 


-9.9 


755 










Phoenix 


dSph 


-9.8 


395 




1 




Saviane et al. 2009 


Sculptor (E351-G30) 


dSph 


-9.8 


87 










Cassiopeia (An VII) 


dSph 


-9.5 


690 










Tucana 


dSph 


-9.6 


870 










Sextans 


dSph 


-9.5 


86 










Carina (E206-G220) 


dSph 


-9.4 


100 










Draco (DDO 208) 


dSph 


-8.6 


79 










Ursa Minor 


dSph 


-8.5 


63 










Canes Vcnatici I* 


dSph 


-7.8 


220 










Leo T* 


dSph 


-7.1 


420 










Ursa Major* 


dSph 


-6.7 


100 










Canis Major Dwarf* 


Irr 




7.6 










Canes Vcnatici II* 


dSph 


-5.8 


150 










Bootes* 


dSph 


-5.8 


60 










Ursa Major II* 


dSph 


-3.8 


30 










LG outskirts 
















GR8 


dSph 


-11.8 


2200 







Magrini et al. 2005 




Antlia 


dSph 


-15.8 


1330 










NGC3109 


dlr 


-15.8 


1330 


18 


20 


Corradi et al. 2006 


Pcna et al. 2007 


Sextans B 


dlr 


-14.3 


1600 


5 


5 


Magrini et al. 2000 




Sextans A 


dlr 


-14.2 


1320 


1 


1 


Magrini et al. 2003 




EGB0427+63 


sir 


-10.9 


2200 











Notes: 

1 This presentation. 2 Hcrndndcz-Martmcz et al. (2009). 



to occur at lower mctallicities. With Type I PNe denned as He/H > 0.10 and N+/0+ > 
0.25, at least 21 PNe from both clouds were found to agree with the N/O requirement 
but low in He (He/H <0.10). These PNe, named as Typelz, were also found by Reid 2007 
in an independent study which used the same method for abundance determination. By 
plotting N/O versus O, Leisy & Dennefeld (2006) found that the destruction of O in high 
mass stars is greater in the LMC and SMC than in our Galaxy. A similar plot for PNe in 
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the SMC by Shaw et al. 2010 gives the impression that there is no substantial decrease 
in N/O with increasing O. These same authors did not find Type I PNe in their observed 
sample but made mention that SMP PN MG8, a non Type I PN, had an estimated cen- 
tral star mass of O.88M implying a main sequence mass > 3.5M Q , which is close to the 
minimum mass where SMC metallicites would allow hot bottom burning to convert C to 
N. 

4.3. Infrared observations 

The growing use of Spitzer SAGE data (Meixner et al. 2006) is enabling far more detailed 
analysis of compact and point-like objects such as PNe in the Magellanic Clouds. With 




Figure 2. Multi- wavelength data for LMC PN SMP47. a-c = HST data, d = UK Schmidt 
Telescope deep ground-based imaging, e = ATCA + Parkes with 8.64 GHz overlaid contours, 
f-k = Spitzer Space Telescope infrared data, L = optical spectroscopy from 2dF on the AAT. 

the increasing use of this data comes a new understanding of the evolutionary properties 
of PNe and the effect that the metallicity of the host system imposes on the observed 
spectra. Using the brightest 233 LMC PNe, Hora et al. (2008) found low detection rates, 
for example, only 119 sources in all 4 IRAC bands, 185 in 2 bands and 28 in only 1 
IRAC band. In MIPS they found 161, 18 and matching the 24, 70 and 160/zm bands 
respectively. The larger number of detections at 24/im probably represents the high 
sensitivity of this MIPS array and the peak of the SED for many PNe which correspond to 
the temperature of ^100K for dust inside the ionised zone, trapped by Lya photons. Low 
detection rates in other bands may indicate that dust and PAHs are evolving according 
to the metallicity, temperature and UV luminosity of the nebula. Compared to Galactic 
PNe adjusted to the LMC distance, the LMC sample has similar 3.6/im to 8/j,m line 
strengths but the sensitivity cutoff for SAGE passes roughly through the middle of the 
Galactic distribution indicating that a large sample of LMC PNe (highly evolved and 
faint ones) are not likely to be detected. Crowding and completeness limits are likely to 
contribute to the non-detections, however using those that could be detected Hora et al. 
(2008) began to find curious differences between Galactic and LMC PNe. 

Areas of IRAC colour-colour plots which are normally associated with galaxies, AGNs, 
and young stellar objects (YSOs) due to star-forming regions in the LMC, are also 
occupied by a large proportion of LMC PNe. In fact, LMC PNe occupy and extend 
all the way to the region populated by stars at 0-0 colour points. Such PNe (eg. J21 and 
J22) are dominated by their central stars indicating that they are highly evolved PNe 
in which the gas is recombining, thereby eliminating other characteristic MIR processes. 
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A MIR survey, isolated from other spectral regimes, may well mistake these objects for 
YSOs, and this is likely to have already occurred. Worse still, at least one attempt to 
classify or re-classify MC PNe using J, Y and K s colours alone has lead to confusion and 
errors. It is clearly risky to attempt to identify MC PNe when IR detections are below 
SAGE limits and optical spectra are not obtained. 

Strong PAH features are typical of young PNe. In more evolved PNe the PAHs are 
dissolving and H recombination lines dominate the near-IR emission (Hora et al. 2008). 

Colour comparisons to 8^im emission show peculiarities of LMC PNe that may be 
linked to low metallicity and/or evolution of the nebulae. Since the 5.8/zm and 8/im 
bands both include PAH emission features, the strengths are well correlated but when 
comparing 4.5/zm to 8/zm for the same objects the greater spread may be due in part to a 
lack of PAHs and much fainter emission from warm dust at 4.5/im. Going further to the 
blue and comparing 3.6^m to 8/xm colours separates PNe such as SMP 11 which occupy 
an area known for extreme AGB stars. While the spectrum of SMP 11 clearly indicates a 
PN, its MIR colours suggest a post- AGB object (Bernard-Salas et al. 2006). Comparing 
8^m versus J-8^m, LMC PNe are appear brighter and redder than the Galactic PNe, 
but this may include a sampling effect since the faintest LMC PNe were not included. 

Great care needs to be taken when interpreting the spectral energy distribution (SEDs) 
for MC PNe as the shape of the SED is strongly influenced by the mass and evolutionary 
state of the PN. When MIR bands for bright, young MC PNe with detectable emission 
at 70/im are plotted, the SEDs all move in a similar line. This is not the case for more 
evolved objects where Hora et al. (2008) show a much more diverse range of SEDs, many 
of which decrease towards 8^m and even decrease from 8^m to 24^im. Similar features 
were found by Stanghcllini et al (2007), who, with the advantage of complete MIR spectra 
found that many LMC PNe show an absence of grain emission and can be categorised as 
"featureless" . PNe with a low dust continuum and the absence of features have not been 
observed in the Galaxy. PNe with a prominent silicon carbide feature at 11/xm are those 
with carbon-rich dust, categorised as type CRD. PNe with crystalline silicate features at 
23.5, 27.8 and 33.5^m have oxygen-rich dust and may categorised as type ORD. 

With these 3 categories identified, the IR luminosity versus carbon abundance shows 
that the carbon-rich PNe are also among the most luminous. Small samples indicate 
they are either round or bipolar core in morphology. Oxygen-rich PNe are also bright 
but may be as much as 1 log C/H lower in carbon abundance than carbon-rich PNe. 
They appear to be bipolar or quadrupolar in morphology. Featureless PNe extend from 
the mid-luminosity region occupied by carbon-rich PNe but then show a steady decline 
in both luminosity and carbon abundance. They appear to be either round or bipolar 
core when carbon-rich but mainly bipolar core when carbon poor. The dust temperature 
also decreases with decreasing IR luminosity where featureless PNe dominate the lower 
end on both scales. These findings show that LMC PNe show extreme changes in their 
MIR and JHK features with evolution. Future work with larger samples will enable clear 
definitions to be formulated. 
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